Gamow- Teller transitions to ^^Cu measured using the ^''Zn(t,^He) reaction 
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The ®''Zn(t,^He) reaction has been studied using a secondary triton beam of 115 MeV/nucleon 
to extract the Gamow- Teller transition-strength distribution to ®^Cu. The results were compared 
with shell-model calculations using the p/-shell effective interactions KB3G and GXPFIA and with 
existing data from the ^'^Zn{d,^}ie) reaction. Whereas the experimental results exhibited good 
consistency, neither of the theoretical predictions managed to reproduce the data. The implications 
for electron-capture rates during late stellar evolution were investigated. The rates based on the 
theoretical strength distributions are lower by factors of 3.5-5 compared to the rates based on 
experimental strength distributions. 

PACS numbers: 21.60.Cs,25.55.-e,25.55.Kr,26.30.-k,26.50.-|-x 



I. INTRODUCTION 

Weak reaction rates play an important role in late stel- 
lar evolution In particular, nuclear electron captures 
(EC) affect the pre-explosion evolution of core-collapse 
supernovae [2|, where the high Fermi energy of degen- 
erate electrons lift Q-value restrictions on EC that are 
present at terrestrial electron densities. The EC reac- 
tions modify properties of the stellar core by neutroniz- 
ing nuclear matter and reducing the electron abundance 
(Ye)- As a consequence, pre-collapse dynamics are mod- 
ified due to the reduced outward pressure. In addition, 
energy and entropy are reduced by neutrino emission as- 
sociated with the weak transitions Fuller, Fowler and 
Newman (FFN) first showed that EC in core-collapse su- 
pernovae are predominantly allowed Gamow- Teller (GT) 
transitions in pf- and sd^-shell nuclei 0, H, 0, 0] ■ They 
provided the first parameterized EC rates for p/-shell 
nuclei, treating nuclear excitations in an independent- 
particle model (IPM) up to A = 60 and assigning a flat 
GT response to heavier nuclei. 

It has since then been understood that residual in- 
teractions between constituent nucleons move the cen- 
troid of the GT strength (B(GT); here defined such that 
B(GT)=3 for the /3-decay of a free neutron) distribu- 
tion and spread it over many more states in the daugh- 
ter than assumed in the IPM used by FFN Shell- 
model calculations employing effective interac- 
tions [ill, [l3| with parameters that are fitted to exper- 
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imental data predict significantly reduced low-lying GT 
strengths and hence reduced EC rates compared to the 
estimates by FFN. Core-collapse models using these re- 
vised rates [1, [l^] indicate that the pre-explosion evolu- 
tion is strongly impacted by the reduced EC rates. It 
is, therefore, crucial that the theoretical models that aim 
to describe GT strength distributions are tested against 
data. The features of the GT strength distributions at 
low excitation energies are particularly important [l[ . At 
any given time prior to the collapse, a varying ensem- 
ble of nuclei plays a significant role [1, Therefore, 
it is important to systematically test theory against ex- 
periment over a broad mass range so that the effective 
interactions can be improved. 

Investigations of charge- exchange (CE) reactions at in- 
cident energies > 100 MeV/nucleon are the pre- 
ferred method for this purpose since, in contrast to /3- 
decay studies, the full excitation-energy region of inter- 
est for the astrophysical applications can be probed. In 
the ATz = +1 direction, which is of direct relevance for 
electron-capture, the {n,p) (see e.g. [13, [13 )j (d,^He) 
(see e.g. [llllO,!!!) and {t,^Re) ]^ probes have been 
employed. Experimental studies of AT^ = —1 probes, 
such as {p,Ti) and (3He,t) have also been used to extract 
GT strength distributions of relevance for EC, under the 
assumption of isospin symmetry [2^, .24] . The foundation 
for the experimental work is the empirical proportional- 
ity between B(GT) and differential cross section for GT 
transitions at vanishing momentum transfer, discussed in 
more detail below. 

In this work, we employed the ^"'Zn(t,3lle) charge- 
exchange reaction at 115 MeV/nuclcon incident triton 
energy to extract the GT transition strength to ^''Cu. 
The experimental results were compared to shell-model 
calculations employing the effective interactions KB3G 
[m and GXPFIA [l|. EC rates on the ^'^Zn ground- 
state were calculated as a function of stellar temperature 
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and density to investigate the impact of the differences 
between the measured and theoretical GT strength dis- 
tributions. We also compare our results with those of 
a recent ^^Zn((i,^He) experiment, which was aimed at 
extracting information on GT strengths in relation to 
double-/? decay transition matrix elements [2ll |. 



II. EXPERIMENT 

A secondary triton beam was produced by fast frag- 
mentation of a 150-MeV/nucleon i^O®"'" primary beam 
on a thick (3500 mg/cm'^) ^Be production target at Na- 
tional Superconducting Cyclotron Laboratory's (NSCL) 
Coupled Cyclotron Facility (CCF) The A1900 

Fragment Separator (26j was set to a magnetic rigid- 
ity Bp = 4.8 Tm, corresponding to 115 MeV/nucleon 
triton energy. A slit at the intermediate image of the 
A1900 was used to limit the momentum acceptance to 
dp/p — ±0.21%, corresponding to ~3 MeV triton kinetic 
energy spread. The isotopic purity of the triton beam was 
85%, with ^He being the principle contaminant, and the 
average intensity at the reaction target was ~ 3 x 10^ 
s"""^. Details of the production of the secondary triton 
beam can be found in Refs. [27l. [28|. 

A 9.84-mg/cm^ thick ^''Zn foil, with an isotopic purity 
of 99.6% was placed at the pivot-point of the S800 spec- 
trograph The analysis beam line to the target was 
operated in momentum dispersion-matching mode [s^l to 
optimize the energy resolution. The S800 spectrograph 
was set at 0° laboratory scattering angle. Although it is 
at present possible to monitor the absolute triton beam 
intensity at the target with good precision (uncertainty 
±5%), at the time of the ®^Zn(i,'^He) experiment sys- 
tematic errors were significantly higher. Therefore, us- 
ing the ^^C(i,'^He) reaction, the excitation of the strong 
GT transition from the i^C 0+ groundstate to the ^^B 
1"^ groundstate with known cross section [13, [3l[ served 
to calibrate the absolute beam intensity. A deuterated 
polyethylene foil (CD2) with a thickness of 9.10 mg/cm^ 
was used for that purpose. Minor variations of the beam 
intensity during ^"'Zn runs due to fluctuations in the ^^0 
primary beam intensity were monitored using a Faraday 
bar placed in the first dipole magnet of the A1900. 

^He ejectiles were momentum analyzed in the S800 fo- 
cal plane using two cathode-readout drift chambers 
(CRDCs) as tracking detectors and two plastic scintilla- 
tors {El and E2) to measure energy loss and event time- 
of- flight (TOF). The El signal was the data-acquisition 
master trigger and TOF start. The CCF radio frequency 
(RF) signal was the TOF stop. Ion energy loss in El 
and event TOF allowed for unambiguous identification of 
•^He events in the S800 focal plane. A 5*'*-order inverse 
ion-optical transfer matrix, calculated based on measured 
magnetic field maps with the code COSY Infinity [33t . 
was used to reconstruct the ^He momentum and scat- 
tering angle (8iab('^He)). The excitation energy in ^''Cu 
(Exi^'^Cu)) was calculated in a missing-mass calculation. 
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FIG. 1: Differential cross section of the '^''Zn(t,^He)'^''Cu CE 
reaction, plotted as a function of excitation energy (Ex) in 
^*Cu. (a) The spectrum obtained for the full solid angle cov- 
ered in the experiment, (b) The spectrum gated on forward 
angles Oiab(^He)< 0.8°. (c) The spectrum gated at backward 
angles 2.4° < eiab(^He)< 3.2° . 



In Fig. [H^a), the excitation energy spectrum up to 15 
MeV for the ^■*Zn(t,^He) reaction is displayed. The exci- 
tation energy resolution was 280 keV (FWHM). Events 
with a scattering angle of up to 4° were included in the 
analysis. The scattering angle resolution was 10 mrad 
(FWHM). 

GT transitions are associated with spin-transfer 
{AS =^ 1) and zero units of angular-momentum trans- 
fer {AL — 0). In the case of the ^^Zn target (ground- 
state = 0+), 1+ states are populated in ^"'Cu. The 
GT differential cross section peaks at 0°, which separates 
them from transitions associated with single or multiple 
units of angular-momentum transfer. Since GT transi- 
tions are in general not separated in excitation energy 
from other transitions, a multipole decomposition analy- 
sis (MDA) is required to isolate the GT contribution to 
the excitation-energy spectrum. In addition, 0+ 1+ 
transitions contain a AL — 2, AS = 1 component. Near 
0°, it is usually small compared to the AL — 0, A^* = 1 
contribution, but the two contributions can only be sepa- 
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rated through a MDA. The MDA performed for this work 
is described in more detail in section Hi Al 

Using the Eikonal approximation, Taddeucci et al. [s^l 
showed that for i^bcam ?J 100 MeV/nucleon, the differen- 
tial cross section at vanishing momentum transfer (q — 0) 
is proportional to B(GT): 



da 

In 



aGTS(GT), 



(1) 



where ctgt is referred to as the unit cross section. This re- 
lationship, first studied for the {p,n) reaction, has since 
been employed for a variety of CE probes. Because of 
the similarity of the ('^He,t) and (tj'^He) probes, the most 
relevant for the present work is the recently empirically 
established [3^1 relationship between the unit cross sec- 
tion and target mass number for the ("^HCjt) reaction at 
140 MeV/nucleon, as discussed in detail in Section HTB] 
The proportionality expressed in Eq. [T]is not perfect. 
The main source of proportionality breaking is the inter- 
ference between AL = 0, AS" = 1 and AL = 2, AS" = 1 
amplitudes due to the non-central tensor component of 
the projectile-target interaction [H, HI]. The effects are 
difficult to quantify experimentally, but can be estimated 
theoretically, as is done in Section Hi CI 
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FIG. 2: Differential cross section from the bin centered at 
3.125 MeV seen in Fig. [Ub). 



A. Multipole Decomposition Analysis (MDA) 

In Fig. [1] the experimental spectra over the full open- 
ing angle covered (Fig. [Ha)) and gated on scattering 
angles at or below 0.8° (Fig. [TJb)) and between 2.4° and 
3.2° (Fig. [IJc)) are shown. The energy bin size is 250 
keV, which is close to the experimental energy resolu- 
tion. For states below the threshold for decay by particle 
emission (the binding energies for protons and neutrons 
are 7.2 MeV and 7.9 MeV, respectively), the line-width 
is much less than the experimental resolution and below 
4° scattering angle most of the yield thus falls within 
one or two energy bins. States associated with AL = 
peak at forward scattering angle and are thus relatively 
strong in Fig. [^b) and weak in Fig. [T]^c). The peak 
seen in the bin at an excitation energy of 3.125 MeV is 
the clearest example. In contrast, the structure seen in 
the spectra at ~ 10 — 12 MeV has a relatively large cross 
section at backward scattering angles (Fig. [He)) but is 
not easily identified at forward angles (Fig. [Hb)). This 
is a characteristic of dipole transitions (AL = 1), which 
peak at 8iab('^He)=2 — 3° and have a local minimum at 
0iab(^He)=O°. Differential cross sections for transitions 
of higher angular-momentum transfers are relatively flat 
at forward angles, with a less-pronounced maximum at 
0iab(^He)=3 — 4°. It is further noted that, apart from 
the GT transitions, the only other significant AL = 
contributions to the excitation energy spectrum in the 
= +1 direction stem from isovector (spin- flip) giant 
monopole resonances. They have negligible cross sections 
below w 10 MeV In this work, the B(GT) in 



^ Cu is extracted only up to 6 MeV. Consequently, the 
assignment of AL = contributions to GT transitions is 
unambiguous. To extract the magnitude of the AL = 
cross section in each 250-keV wide bin, a MDA was per- 
formed. In Fig. [21 an example of the procedure is shown 
for the bin centered at 3.125 MeV. The '^He angular dis- 
tributions were generated in the Distorted Wave Born 
Aroroximation (DWBA) with the DWBA code FOLD 
[3^1; transitions involving AL = 0, AL = 1 and AL — 2 
to ^''Cu were considered. Transitions with AL > 2 are 
suppressed at these beam energies and their angular dis- 
tributions below 4° scattering angle are quite similar to 
the L = 1,2 shapes already used in the MDA. The ef- 
fective NN interaction by Love and Franey at 140 MeV 
[39I . l40j was double- folded over the transition densities of 
the t-^He and ^^Zn-^'*Cu systems to generate the form 
factor. A short-range approximation [3^ was used to 
describe exchange processes. Binding energies of single- 
particle wave functions in ®^Zn and ^^Cu were calculated 
in the code OXBASH t41|] using the Skyrme SK20 in- 
teraction [i^l. The ^He and densities were obtained 
from Variational Monte-Carlo results [i^. One-body 
transition densities for GT transitions were generated in 
NuShellX [44j using the GXPFIA [12] interaction. 

The shape of the AL = angular distributions do not 
vary significantly for transitions to different states. For 
transitions associated with AL = 1 or 2, one-body transi- 
tions densities were generated in a normal-mode formal- 
ism using the code NORMOD [4^. Optical-potential 
parameters for both entrance and exit channels were 
taken or deduced from ^He elastic-scattering data on ^*Ni 
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p^ . The ^He-^^Cu optical-potential parameters were - 
35.16 MeV (-44.43 MeV), 1.320 fm (1.021 fm), and 0.840 
fm (1.018 fm) for real (imaginary) Woods-Saxon well- 
depth, radius and diffuseness, respectively. For the t- 
^^Zn system, the above well-depths were scaled by a fac- 
tor 0.85, following the procedure in Rcf. 47||. 

Pair-wise combinations of calculated angular dis- 
tributions associated with different units of angular- 
momentum transfer were used to fit the experimental 
angular distribution for each excitation-energy bin. Only 
the normalizations of the calculated angular distributions 
were taken as free fit parameters. Due to the limited an- 
gular coverage, fits with more than two components did 
not improve the reduced value of the fit, nor change 
the cross section of the AL = contribution by more 
than a fraction of the statistical errors. The fit with 
the combination of angular distributions minimizing the 
reduced value was chosen as the best description of 
the data. The angular distribution of the bin centered 
around 3.125 MeV shown in Fig. [2] was best described 
by a strong J'^ = 1+ component in combination with a 
weaker = 2+ contribution. Statistically significant 
AL = contributions to the spectrum were found up 
to an excitation energy of 6 MeV. In each energy bin, 
the fitted AL = curve was used to determine the cross 
section at 0°. The fit error was used to determine the un- 
certainty in the extracted 0° cross section. In some bins 
no significant AL = contribution was found. Although 
it cannot be excluded that small AL — contributions 
are present in these bins due to the limited statistics, we 
report no Gamow- Teller strength for such cases. 

The 0° cross section for a transition at Q < must 
be extrapolated to Q = to reach the zero linear- 
momentum transfer (g = 0) limit required for application 
of Eq. [T] The extrapolation factor was determined using 
the relation: 



da 



9^0 



^(Q.o°) 



DWBA 



da 

dh 



(Q,0°) 



(2) 
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where the subscript "DWBA" denotes calculated cross 
sections in DWBA and the subscript "exp" denotes the 
measured cross section. The DWBA calculations were 
performed for 100 reaction Q-values so that the ratio of 
DWBA cross sections in Eq. [2] could be established as a 
smooth function of excitation energy. 



B. Application of the (''He,t) Unit Cross Section 

From the extracted cross sections at q = 0, abso- 
lute B(GT)s were determined by dividing the measured 
cross section by the unit cross section acx- The unit 
cross section is often calibrated directly using a transi- 
tion for which the B(GT) is known from a /3-decay mea- 
surement. In principle, this is possible for the case of 
^■^Zn^^^Cu as well, since the log/t for the groundstate- 
to-groundstate transition is known 48]. The measured 



value is log/t=5.301±0.006, which, following the formal- 
ism of Ref. corresponds to a B(GT) for the /3-decay 
of ^"^Cu of 0.0192 ± 0.0004. Taking into account that 
B{GT)a^t - {2Jb + l)/(2Ja + l)B{GT)b^a, the B(GT) 
for the 64Zn-^64Q^ transition is 0.058±0.001. This direct 
calibration of the unit cross section was used in the anal- 
ysis of the ^^Zn((i,^IIe) experiment [21*1. However, due to 
the limited energy resolution and statistics in the present 
(i,^He) data set, the relatively weak groundstate transi- 
tion could not be resolved from the 1+ excited state at 
344 keV Moreover, as discussed in the next section, 
the calibration of the unit cross section using a /3-decay 
strength to a weak state can lead to significant errors if 
the proportionality breaking is large for the particular 
transition used due to effects of the tensor interaction 
(see e.g. the case ^^Ni 

As an alternative method, we use the empirically es- 
tablished mass-dependent trend line for the unit cross 
section of the (^He,i) probe [35j: 



CTGT = 109 A 



-0.65 



(3) 



where A is the target mass number. It was recently found 
that the unit cross sections for the {t,^lle) reaction in 
the case of A=12,13 and 26 were consistent within error 
margins to that of the (■^He,t) reaction 31], in spite of the 
minor difference in beam energy (115 MeV/nucleon for 
(t,3He) versus 140 MeV/nucleon for (^He,^)). Although 
future (tj^He) experiments over a wider mass range are 
important to further establish the consistency of the unit 
cross sections between the two inverse reactions, it seems 
likely that the close correspondence will be maintained 
for higher mass numbers. 

In summary, Eqs. [1] [21 and [3] were used to calculate 
the absolute B(GT) of transitions to states in ^''Cu based 
on the extracted AL = partial cross sections. Fig. [3] 
shows the resulting B(GT) distribution in ^''Cu. The 
transition to the 3.125 MeV bin was identified as having 
the highest B(GT). The B(GT) values in the first two 
250-keV wide bins were adjusted to account for the fact 
that the GT strength for the groundstate transition is 
known. Hence, the B(GT) value in the first bin was ad- 
justed to the value obtained from /3-decay. The remaining 
strength observed experimentally in that first excitation- 
energy bin was shifted to the second excitation-energy 
bin under the assumption that, due to the limited reso- 
lution, some of the higher-lying strength was detected at 
too low an excitation energy. 

The errors presented in Fig. [3] are dominated by the 
statistical error, with a small contribution from uncer- 
tainties due to the assumptions used in the fitting pro- 
cedure. There is an uncertainty in the overall scal- 
ing of about 10% due the relative normalization to the 
^^C(i,^He)^^B(l+) reaction and the application of the 
unit cross section as calculated in Eq. [3] 
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FIG. 3: The GT strength in ''''Cu, extracted up to 6 MeV in 
excitation energy. 



C. Analysis of Systematic Errors 

As mentioned above, the proportionality between 
B(GT) and the cross section at g = is not perfect, due 
to interference between AL = 0, AS* — 1 and AL = 2, 
AS* = 1 ampHtudes. The latter amplitude is mediated 
through the tensor component of the effective NN inter- 
action. In general, it is difficult to evaluate these er- 
rors on a transition-by-transition basis, since it requires 
an accurate knowledge of the contributions from individ- 
ual Ip — l/i components to the transition density [23 |. 
However, as shown in Ref. [s^, a theoretical study of 
the proportionality breaking for many shell-model states 
can reveal the associated uncertainties as a function of 
B(GT) and the same procedure was carried out for the 
^^Zn(t,'^He) reaction and is described in this section. 

To evaluate the systematic error due to the interfer- 
ence, the one-body transition densities of 100 GT shell- 
model transitions generated with the GXPFlA [12] in- 
teraction (all states predicted below about 6 MeV) were 
used to calculate differential cross section in DWBA. 
From there on, the theoretical cross sections were treated 
as if they were data and B(GT)s were extracted as de- 
scribed above. Since the extraction assumes a common 
unit cross section for all transitions, the deviation be- 
tween the extracted strengths (i3(GT)DWBA) and the 
values calculated in the shell model (i3(GT)sM) provides 
a measure for the amount of proportionality breaking 
due to the interference and other effects included in the 
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FIG. 4: A theoretical estimation of the relative systematic 
error in the B(GT+) extraction using Eq. [1] plotted as a 
function of the shell-model strength. 



DWBA calculations. A relative systematic error was de- 
fined as: 



Rel.sys.err. 



S(GT)dwba - S(GT) 



SM 



S(GT)s 



M 



(4) 



Fig. 2] shows the relative systematic error for the 100 
^^Zn(i,'^He)^'*Cu(l"'") transitions generated, plotted as a 
function of the B(GT) determined in the shell-model. On 
average, the relative systematic error increases with de- 
creasing B(GT). An approximate relationship between 
the width of the error distribution and the B(GT) was 
extracted from these points: 



'^rel.sys.err. ~ 0-03 - 0.033 x ln(B(GT)). 



(5) 



It is noted that this relationship is very similar to that 
found for the case of ^^Mg('^He,t) in spite of the fact 
that the shell-model spaces and interactions employed 
are very different. 

The procedure for estimating the breaking of the unit 
cross section in the ^^Zn case was performed a second 
time after switching off the tensor contributions to the 
effective NN interaction [s^, |4^. The relative system- 
atic error as defined in Eq. 3] became much smaller; 
(Tj.gj gj.j. < 5%, even for the weakest transitions. 
These remaining errors are due to effects such as ex- 
change contributions and small differences in the shape 
of the transition densities for different GT transitions, 
but are much smaller than the effects due to the tensor 
interaction. 

There are outliers in the relative systematic error dis- 
tribution. For example, for the state seen in Fig. 0] at 
i3(GT)sM ~ 0.03, which corresponds to the first excited 
state in the theoretical spectrum, the proportionality is 
broken by « 50%. If such a strong breaking were to hap- 
pen for the groundstate transition and that transition 
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FIG. 5: (color online) Running sum of B(GT) in ''"Cu, plotted 
as a function of excitation energy, for B(GT) extracted with 
the (t,^He) reaction (black, thick cross hatches) and with the 
(dj^iie) reaction (red, thin cross hatches) [2ll |. 



were to be used for the overall calibration of the unit 
cross section, instead of the present procedure, the same 
large error would be introduced for all transitions. Ap- 
plying the (^He,t) unit cross section as described in Eq. 
[3] to (t,^lle) cross section data to extract B(GT) guards 
one against this potentially large systematic error. 



D. Consistency with the (d,^He) Reaction 

In Fig. [5] the running sum of B(GT) as a function 
of the excitation energy in ^*Cu is compared with that 
of the ^'*Zn((i,^IIe) measurement of Grewe et al. [2l| . 
The width of the bands represent the cumulative errors 
in the summations. Note that in Ref. [2l[, B(GT)s 
were extracted only up to 5 MeV. The main features 
of the two distributions are consistent. A difference of 
K, 10% in the overall strength up to an excitation en- 
ergy of 5 MeV is found, but the deviations are within 
error margins, which, moreover, do not include errors in 
the overall B(GT) normalization for either data set. The 
summed B(GT)s up to an excitation energy of 5 MeV 
are ^-B(GT) = 1.78 ± 0.13 in the present work and 
X; B(GT) = 1.61 ± 0.05 for the results from the {d,'^m) 
experiment. 

The reactions and methods employed to extract the 
B(GT) distribution were quite different for the present 
work and that of Grewe et al. [2l|. Here, the assumption 
was made that the unit cross section from (•^He,t) could 



be employed and the analysis of the (d,^He) data relied 
on the fact that the proportionality between differential 
cross section at g = and B(GT) was not strongly bro- 
ken for the groundstate transition that was used to cali- 
brate the unit cross section. The fact that the results are 
consistent gives confidence that both assumptions were 
reasonable, on the level of « 10%. 



III. COMPARISON WITH SHELL-MODEL 
CALCULATIONS 

The goal of the present work is to test the theoretical 
predictions for weak-transition strengths that are used 
to generate electron-capture rates of relevance for stellar 
evolution. In this section, we compare theoretical and 
experimental GT strength distributions. In the next sec- 
tion, the implications for the electron-capture rates under 
stellar conditions are explored. 

The shell-model has provided the best description of 
low-lying GT strength distributions [o', [l3|. Here, our ex- 
perimental results are compared with shell-model calcu- 
lations using the pf- shell effective interactions GXPFIA 
[l5| and KB3G llf. The parameters of the GXPFIA in- 
teraction have been fitted to reproduce the experimental 
excitation energies and masses (but not GT strengths) 
over a wide range of p/-shell nuclei including A=64. 
KB3G is a slightly modified version of the KBF inter- 
action Their parameters have been deduced from 
experimental data in the lower p/-shell only. The KBF 
interaction has been used to generate weak transition rate 
tables [50] that are commonly used in astrophysical sim- 
ulations. 

The shell-model calculations were performed using the 
code NuShellX [ill , which was slightly modified to pro- 
vide a large number of states of a given spin. Figure 
[6] shows the running sum of B(GT) in ^"'Cu up to an 
excitation energy of 7.5 MeV. Results from the present 
experimental work and from the two theoretical calcula- 
tions using the GXPFIA [l2| and KB3G [III interactions 
are shown. The theoretical calculations have been scaled 
by a p/-shell specific quenching factor of (0.74)^ [5l| to 
account for effects not included in the model space. 

Neither calculation reproduces the experimental 
strength distribution. The distribution for GXPFIA is 
closer to the data, but pushes the strength up too high 
in excitation energy. The same, but more dramatically, 
happens for the calculation using KB3G, although the 
strength integrated up to 7.5 MeV reproduces the ex- 
perimental value quite well. The summed B(GT) up to 
E.^ = 7.5 MeV (a total of 48 states) for the KB3G inter- 
action is E-S(GT)kb3G = 2.02 (a further 10% of that 
value is located at energies up to 10.3 MeV) compared to 
the experimental value of 1.95±0.14 up to that excitation 
energy. The summed strength up to = 7.5 MeV with 
the GXPFIA interaction is E^(GT)gxpfia = 2.65. A 
further 8% of that value is located at higher excitation 
energies, fragmented over many weak states. 
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FIG. 6: (color online) Running sum of B(GT) in '^"Cu, plot- 
ted as a function of excitation energy, as measured with 
the (t,^He) reaction (black, cross hatched), and calculated 
with the shell-model using GXPFIA 12] (red, solid line) and 
KB3G yjl] (green, dashed line) effective interactions. 
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FIG. 7: (color online) Electron-capture rate on ®'*Zn, plotted 
as a function of stellar temperature, at an electron density 
of 10^ g/cm^. Rates based on GT strength from the (f,'^He) 
data (black, solid line), the (d,^IIe) data (blue, dotted line), 
the shell-model calculations using GXPFIA (red, dashed line) 
and KB3G (green, dashed-dotted line) are shown. The yellow 
cross-hatched region marks the relevant temperature range at 
this density in the pre-coUapse trajectory for a ISMgoiar mass 
star. 



IV. APPLICATION TO ELECTRON CAPTURE 
RATES 



A possible reason for the discrepancy between the data 
and the theory is a contribution from the gg/2 orbit. In 
Ref. ^5^], the observation of a collective rotational band 
in ^*Zn, extending from spins 12fi. to 24?i with the band 
head at ~ 8 MeV was reported. That observation was 
explained by a model in which a hole was created in the 
proton /7/2 orbit associated with triaxial deformation. 
One proton and two neutrons then resided in the 179/2 
orbit [H,!!^. Hence, at higher excitation energies one 
can expect contributions from the 179/2 orbit, which are 
not included in the shell-model calculations performed for 
this work to calculate the GT strength distributions. At 
the excitation energies of relevance here, the influence of 
the 39/2 orbit is not so clear. Theoretical studies of shape 
coexistence [sJl indicate that the effects of the 39/2 are 
small in Zn nuclei and in a recent study of the structure 
of ^''Ge [5^, E2 transition rates were well-described us- 
ing the GXPFIA interaction. In summary, the intrusion 
of the (79/2 shell likely plays a role, but not necessarily 
at _Ea;(^'*)Cu< 6 MeV. The strength distribution is also 
affected by the differences in deformation of ^''Zn and 
^'^Cu, which complicates matters further. 



Calculations of the stellar EC rates onto the ^ Zn 
groundstates were carried out following the procedures of 
Refs. using a code by Gupta et al. [5^^. The 

capture rate was calculated over a wide grid of tempera- 
tures (0.01x10^ <T < lOOxlO^K) and electron densities 
(10^ < pYc < lO^^g/cm^) which should cover nearly all 
scenarios relevant to nuclear astrophysics. Capture-rate 
calculations presented here are for two of the above elec- 
tron densities; at pY^ = lO^g/cm'^ and pYc = lO^g/cm^, 
chosen for their relevance along the core-collapse super- 
nova trajectory. Electron densities from pY^ ~ lO^g/cm'^ 
to pYc = lO^g/cm'^ correspond to the conditions present 
in the core from the silicon-burning phase of the pre- 
collapse progenitor star to the conditions present at 
about 0.5 seconds before core-bounce, respectively. This 
time interval during the pre-bounce evolution is where 
electron-captures on nuclei have potentially the largest 
influence on the post-bounce trajectory [sl. |57|. 

Fig. [7] shows the capture rate as a function of stellar 
temperature at pYc — 10''g/cm^. The cross-hatched re- 
gion indicates the interval of relevant temperatures for 
a 15 solar mass, main-sequence progenitor near the sili- 
con burning phase. Capture-rate curves are plotted for 
the B(GT) extracted from the present experiment, that 
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the KB3G rate by a factor of 5 compared to the rate 
based on the present experimental data. This is in spite 
of the improved performance at increased Fermi energy 
(4.7 MeV at pY^ — lO^g/cm^) and higher excitation ener- 
gies demonstrated by the GXPFIA calculation in Fig. [5l 
The reason for this, as pointed out by FFN [1, HI, is that 
while the branching ratio for capture to a given final state 
is proportional to its B(GT), the phase space available 
for electrons with an energy sufficient to excite levels at 
higher excitation energy in the daughter nucleus shrinks 
much more rapidly. Therefore, the total capture rate at 
the densities and temperatures studied here depends sen- 
sitively on the GT distribution at low excitation energies 
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FIG. 8: (color online) Electron-capture rate on ^''Zn, plotted 
as a function of stellar temperature, at an electron density 
of 10^ g/cm^. Rates based on GT strength from the (f,^He) 
data (black, solid line), the (d,^He) data (blue, dotted line), 
the shell-model calculations using GXPFIA (red, dashed line) 
and KB3G (green, dashed-dotted line) are shown. The yellow 
cross-hatched region marks the relevant temperature range at 
this density in the pre-coUapse trajectory for a ISMgoiar mass 
star. 



resulting from the B(GT) values published by Grewe et 
al. [2l|, and those resulting from B(GT) calculated using 
GXPFIA and KB3G shell-model interactions. 

The electron Fermi energies for the density of pY^ — 
lO^g/cm"^ is 0.7 MeV and the EC rates are thus very 
sensitive to the GT strength distribution at low excita- 
tion energies at this density. The capture rate based on 
the two measured B(GT) distributions differs by about 
10% in the window of relevant temperature. This is due 
to slight differences in the experimentally extracted GT 
strength distributions and the small discrepancy in the 
overall normalization discussed above. The level of dis- 
agreement between the rates based on the (tj'^He) and 
((i,^He) data thus serves as a reasonable estimate of the 
systematic error in the capture rate determination from 
experiment. The comparative lack of strength in the 
first MeV of excitation energy predicted by both shell 
model calculations results in a severe underestimation of 
the capture rate. In Fig. [71 the capture rate using the 
GXPFIA interaction is lower by a factor of 3.5 and the 
KB3G result is lower by a factor of 5, relative to the rate 
based on the strengths extracted from the ^'*Zn(t,''He) 
experiment. 

These differences between experiment and theory per- 
sist at higher densities, as shown in Fig. [51 There, the 
two measured capture rates still differ by about 10%; 
the GXPFIA rate remains lower by a factor of 3.5, and 



CONCLUSIONS 



In summary, the GT strengths for transitions from 
^^Zn to ^''Cu have been extracted using the ^'*Zn(t,^He) 
CE reaction at 115 MeV/nucleon. The extraction pro- 
cedure was based on the empirical proportionality be- 
tween the cross section in the limit of vanishing momen- 
tum transfer and B(GT) derived for the ("^HCji) reaction. 
Systematic errors in the extraction due to action of the 
tensor operator in the effective NN interaction have been 
investigated. It breaks the proportionality by 10-20% 
and is strongest for weaker GT transitions. 

In spite of the different method employed to extract the 
B(GT), the experimentally extracted GT-strength distri- 
bution was, within error margins, consistent with that 
obtained from a ^''Zn(fi,^He) experiment [2l|. The ex- 
perimental GT strength distribution has also been com- 
pared to shell-model calculations using the NuShellX 
code employing the GXPFIA and KB3G effective inter- 
actions. Both sets of calculations fail to reproduce the 
data, although GXPFIA describes the strength distribu- 
tion somewhat better than KB3G. 

Since the shell-model calculations underestimate the 
GT strength at low excitation energies, deduced electron- 
capture rates on the ^^Zn groundstate from the theoreti- 
cal strength distributions using the GXPFIA (KB3G) in- 
teraction are too small compared to those deduced from 
the experimental distribution by a factor of 3.5 (5). The 
results indicate that further refinement of the interactions 
are required. It can furthermore be concluded that addi- 
tional data is required to test the theoretical calculations 
and to check if the deviations are systematic in the up- 
per p/-shell. An underestimation of the electron-capture 
rate on a single nucleus is not likely to strongly affect the 
impact of electron capture on astrophysical scenario, but 
systematic underestimation of electron-capture rates in 
a region of the nuclear chart can have significant conse- 
quences. 
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